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ABSTRACT 

In  this  study  the  subtropical  cells  (STC)  in  the  Pacific  Ocean  are  analyzed  using  an  eddy-resolving  ocean 
general  circulation  model  driven  by  atmospheric  forcing  for  the  years  1992-2003.  in  particular,  the  authors 
seek  to  identify  decadal  changes  in  the  STCs  in  the  model  and  to  compare  them  with  observations  in  order 
to  understand  the  consequences  of  such  changes  for  the  equatorial  ocean  heat  and  mass  budgets.  The 
simulation  shows  a  trend  toward  increasing  pyenoclme  volume  transport  at  9°N  and  9  S  across  Ihe  basin 
from  1992  lo  2003.  This  increase  |4.9  ±  L0  Sv  (Sv  -  10h  mV  is  in  qualitative  agreement  with  obser¬ 
vations  and  is  attributed  primarily  to  changes  in  the  interior  ocean  transport,  which  are  partially  compen¬ 
sated  hv  opposing  western  boundary  transports.  The  subtropical  meridional  volume  transport  convergence 
anomalies  in  Ihe  model  pyenoclme  are  found  to  be  consistent  with  anomalous  volume  transports  in  both  the 
observed  and  modeled  Equatorial  Undercurrent,  as  well  as  with  Ihe  magnitude  of  simulated  anomalous 
up  welling  transport  at  the  base  of  the  mixed  layer  in  the  eastern  Pacific,  As  a  result  of  the  increased 
circulation  intensity,  heat  transport  divergence  through  the  lateral  boundaries  of  the  tropical  control  volume 
(defined  as  the  region  between  9°N  and  9°S,  and  from  the  surface  to  <rrt  =  25*3  isopycnal)  increases,  leading 
to  a  cooling  of  the  tropical  upper  ocean  despite  the  fact  that  net  surface  heal  flux  into  the  control  volume 
has  increased  in  the  same  time.  As  such,  these  results  suggest  that  wind-driven  changes  in  ocean  transports 
associated  with  the  subtropical  cells  play  a  central  role  in  regulating  tropical  Pacific  climate  variability  on 
decadal  time  scales* 


I.  Introduction 

The  subtropical  cells  (STCs)  in  the  Pacific  Ocean 
provide  a  mechanism  for  extra  tropica  l  and  tropical 
communication.  This  circulation  is  important  for  sett  ing 
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up  the  Ihertnocline  mean  structure  in  the  Tropics  by 
subducting  water  in  the  subtropics  and  ventilating  ihe 
equatorial  ocean.  As  such,  if  can  play  a  role  in  regulat¬ 
ing  tropical  interannual  climate  variability  since  phe¬ 
nomena  like  El  Nino  are  critically  influenced  by  the 
structure  of  the  mean  background  condition  (e.g,*  Fe¬ 
dorov  and  Philander  2001),  Recently  variations  in  the 
STC  have  also  been  proposed  as  a  possible  mechanism 
for  low-frequency  climate  variability  (longer  than  inter¬ 
annual)  in  the  Tropics,  Gu  and  Philander  (1997)  sug¬ 
gest  that  decadal  variations  in  equatorial  sea  surface 
temperature  (SST)  are  indue  need  by  the  advection  of 
water  mass  anomalies  from  the  subiropics.  Later  stud- 
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iesshow  that  it  is  not  evident  that  water  mass  anomalies 
generated  at  subtropical  latitudes  can  survive  with  any 
appreciable  amplitude  by  the  time  they  reach  equato¬ 
rial  latitudes  (Deser  et  ah  14%:  Schneider  el  aL  1449). 
An  alternative  role  for  the  STCs  in  decadal  climate 
variability  is  that  variations  in  pycnodine  transport  act 
on  the  mean  hydrographic  structure  and  are  primarily 
responsible  for  changes  in  equatorial  SSTs  (Kleeman  et 
ah  1999:  Nonaka  el  al.  2002).  This  latter  hypothesis  is 
supported  by  observational  studies  (e.g..  McPhaden 
and  Zhang  2002,  2004)  and  coarse-resolution  ocean 
general  circulation  model  simulations  (e,g„  Lee  and  Fu- 
kumori  2003:  Capotondi  et  ah  2005). 

The  strength  of  the  STCs  can  be  represented  by  the 
transport  in  its  lower  limb,  namely,  the  volume  trans¬ 
port  of  water  in  the  ocean  pycnodine  layers.  In  obser¬ 
vational  studies,  ocean  subsurface  hydrographic  data 
are  used  to  estimate  the  geostrophic  transport  beneath 
the  surface  layer,  but  data  are  not  adequate  to  define 
decadal  variations  in  western  boundary  current  trans¬ 
ports  that  represent  a  significant  fraction  of  STC  pyc¬ 
nodine  transport  (Liu  1994).  Modeling  studies  (Lee 
and  Fukumori  2003:  Capotondi  et  ah  2005)  also  suggest 
a  partial  compensation  between  the  interior  ocean  and 
western  boundary  pycnodine  transport  variability.  Yet, 
the  representation  of  western  boundary  transport  in 
those  models  is  questionable  because  of  their  coarse 
resolution. 

In  this  study,  we  investigate  variations  of  the  Pacific 
STC  transport  for  the  period  ] 442-2003  in  an  eddy¬ 
resolving  ocean  general  drculation  model.  The  model  was 
forced  with  6-hourly  European  Centre  for  Medium- 
Range  Weather  Forecasts  (ECMWF)  reanalysis  and 
operational  analysis  surface  fluxes.  Our  study  has  three 
specific  purposes.  First,  the  observational  study  by 
McPhaden  and  Zhang  (2004)  shows  an  upward  trend  in 
the  STC  transport  in  the  Pacific  from  the  early  1990s  to 
early  2000s,  indicating  that  this  circulation  has  re¬ 
bounded  from  a  less  active  stale  that  has  lasted  since 
the  1976/77  regime  shift  in  the  Pacific  basin.  As  ac¬ 
knowledged  by  McPhaden  and  Zhang  (2004),  there  are 
considerable  uncertainties  in  their  calculations  because 
the  observations  are  sparse  in  space  and  time,  the  geo¬ 
strophic  method  requires  assumptions  about  a  level  of 
no  motion,  and  decadal  variations  in  western  boundary 
current  transport  cannot  be  reliably  estimated.  A 
model  simulation  where  all  ocean  fields  are  available  at 
high  resolution  provides  a  valuable  dataset  for  improv¬ 
ing  our  understanding  of  the  relevant  processes. 

Second,  we  seek  to  understand  the  relationships  be¬ 
tween  the  meridional  transport  of  the  STCs  in  the 
pycnodine,  the  Equatorial  Undercurrent  (EUC),  and 
the  upwelling  into  the  mixed  layer  in  the  eastern  tropi¬ 


cal  Pacific.  In  the  context  of  ventilated  thermoelme 
theory  (Luyten  et  al.  1983).  all  of  these  transports  can 
be  viewed  as  branches  of  the  STC:  water  subducted  in 
the  subtropical  outcrop  region  eventually  joins  the 
EUC  after  an  intricate  pathway  through  the  ocean  in¬ 
terior  and  the  western  boundary  (Pedlosky  1987,  1988: 
Liu  et  al  1994).  The  water  is  then  transported  east¬ 
ward  by  the  EUC,  and  finally  upwells  to  the  surface 
layers  in  the  eastern  basin.  However,  very  few  studies 
have  addressed  the  issue  of  decadal  variations  in  EUC 
transports  and  how  they  relate  to  STC  variability 
(Chang  et  al.  2001).  Here  we  use  long  lime  series  data 
from  the  Tropical  Almosphere-Oeean/Triangle 
Trans-Ocean  Buoy  Network  (TAO/TRITON)  array 
(McPhaden  et  ah  1948)  to  validate  model  EUC  varia¬ 
tions  and  relate  EUC  transports  to  other  branches  of 
the  Pacific  STCs. 

The  third  purpose  is  to  examine  how  the  upper-ocean 
heat  balance  in  the  tropical  Pacific  is  influenced  by  the 
volume  transport  variability  of  STC  on  decadal  time 
scales.  Previous  studies  (e.g.,  Hazeleger  et  ah  2004) 
have  suggested  that  total  poleward  heal  transport  in  the 
tropical  ocean  is  controlled  by  local  wind-driven  over¬ 
turning  and  gyre  circulations,  as  well  as  by  extra  tropica  I 
buoyancy-flux-driven  overturning  that  reaches  into  the 
Tropics.  In  this  study,  we  will  examine  the  relationship 
between  anomalous  pycnodine  volume  transport,  its 
associated  heat  transport,  and  SST  using  an  idealized 
model.  We  will  also  examine  the  relative  contributions 
of  circulation  and  temperature  anomalies  to  the  heat 
transport.  Since  modeling  assessments  of  role  of  ocean 
circulation  on  the  upper-ocean  heat  budget  depend  on 
model  physics,  particularly  mixing  para  mete  rizati  on 
our  results  will  be  compared  with  available  modeling 
studies  to  indicate  model  sensitivity. 

The  paper  is  organized  as  follows:  In  section  2  we 
introduce  the  ocean  general  circulation  model  and  the 
atmospheric  forcing  data  used  in  this  study.  The  results 
in  terms  of  STC  volume  transports,  their  decadal  vari¬ 
ability  and  relationships  to  the  EUC  and  equatorial  up- 
welling,  and  impacts  on  upper-ocean  heat  balance  in 
the  Tropics  are  presented  in  section  3.  In  section  4  we 
compare  some  of  our  results  to  other  recently  published 
modeling  results  and  discuss  this  particular  model's 
limitations.  We  summarize  our  main  conclusions  in  sec¬ 
tion  5, 

2.  Model 

The  model  used  in  this  study  is  the  Hybrid  Coordi¬ 
nate  Ocean  Model  (HYCGM)  (Bleck  2(X)2:  see  infor¬ 
mation  online  at  http://hycom.rsma$.miami,edu).  It  cov¬ 
ers  the  Pacific  Ocean  from  20°S  to  66°N.  Using  Merca- 
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tor  projection,  the  horizontal  grid  resolution  is  0.08°  in 
longitude  by  0.08Dcos(<i> )  in  latitude,  where  <5  is  the 
latitude.  This  gives  a  midlatitude  grid  spacing  of  ap- 
proximate  ly  6  km  and  a  tropical  spacing  of  just  under  9 
km.  It  has  20  vertical  coordinate  surfaces  that  are  typi¬ 
cally  isopycnal  in  the  open,  stratified  ocean  buL  makes 
a  dynamically  smooth  transition  to  terrain-following  co¬ 
ordinates  in  shallow  coastal  regions  and  to  c-level  co¬ 
ordinates  in  the  mixed  layer  and/or  unstratified  waters. 
The  hybrid  coordinate  is  obtained  via  a  fully  general 
continuity  equation  that  allows  an  arbitrary  partitioning 
between  density  coordinates  and  depth  coordinates  on 
a  time  step  by  time  step  basis.  Specific  details  about  the 
numerics  can  be  found  in  Bleck  et  ah  (2002),  and  ad¬ 
ditional  case  studies  can  be  found  in  Chassignet  et  al. 
(2003). 

HYCOM  is  configured  with  several  different  mixed 
layer  formulations  (see  Halliwell  2004,  for  more  de¬ 
tails),  but  the  A -profile  parameterization  (Large  et  al, 
1994)  is  used  exclusively  here.  The  bathymetry  is  based 
on  a  merging  of  the  global  5'  Digital  Elevation  Model 
(ETOP05)  (National  Oceanic  and  Atmospheric  Ad¬ 
ministration  1986)  and  Smith  and  Sandwell  (1997)  data 
with  numerous  hand  edits.  Closed  boundaries  are  used 
for  the  Indonesian  Through  flow  region  and  in  the  Ber¬ 
ing  Strait.  Along  the  southern  boundary,  the  ocean 
temperature  and  salinity  are  relaxed  to  3D  monthly 
climatological  data  from  the  Generalized  Digital  Envi¬ 
ronmental  Model  (GDEM).  version  3.  (Carnes  2002)  in 
a  3°~wide  buffer  zone.  The  deformation -dependent 
Laplacian  viscosity  factor  is  0,1  and  a  combination  of 
Laplacian  (0.005)  and  hi harmonic  (0,01)  momentum  dis¬ 
sipation  is  used  for  the  large  and  small  scales,  respectively. 

The  model  was  initialized  from  a  0,16°  version  of 
Pacific  HYCOM  and  spun  up  at  0.08°  to  near  equilib¬ 
rium  using  a  climatology  of  winds  and  thermal  for¬ 
cing  based  on  the  1 979-93  ECMWF  Reanalysis 
atmospheric  products.  After  spin  up,  the  simulation  was 
extended  using  the  6-hourly  1979-93  ECMWF 
Reanalysis  and  then  1994— 2003  ECMWF  operational 
forcing.  No  significant  basinwide  changes  in  forcing 
were  detected  across  the  1993-94  time  boundary.  The 
10-m  wind  velocity  components  were  converted  to 
stresses  using  a  stability-dependent  drag  coefficient 
from  Kara  el  al  (2002).  Thermal  forcing  includes  air 
temperature,  specific  humidity,  and  radiative  (short¬ 
wave  and  longwave)  fluxes.  Precipitation  is  also  in¬ 
cluded  as  a  surface  forcing.  Surface  latent  and  sensible 
heat  fluxes  are  calculated  using  "bulk  formulae”  given 
the  above  atmospheric  conditions  and  model  SST, 
Model  sea  surface  salinity  (SSS)  is  relaxed  to  monthly 
GDEM3  SSS,  but  no  SST  relaxation  is  used  by 
HYCOM.  Monthly  river  runoff  is  parameterized  as  a 


Flo.  1.  Meridional  overturning  circulation  in  I  he  lop  300  m  of 
the  model  (Sv),  Solid  (dashed)  lines  indicate  clockwise  (counter¬ 
clockwise)  circulations.  Dashed-dotted  lines  indicate  the  22.3  and 
25.3  trH  surfaces.  Mean  state  for  (a)  Period  1  (defined  in  the  text), 
(b)  Period  2  (defined  in  the  text),  and  (c)  the  difference  bciwcen 
Periods  1  and  2,  (b)  minus  (a),  are  shown. 

surface  precipitation  flux  al  a  select  number  of  loca¬ 
tions.  and  oceanic  turbidity  is  also  taken  into  account 
after  the  work  of  Kara  cl  al.  (2005). 

All  model  output  used  in  this  study  consists  of 
monthly  averages  from  year  1992  to  2003,  T  he  reason 
we  focus  on  this  period  is  to  facilitate  a  direel  compari¬ 
son  with  the  observational  study  of  McPhaden  and 
Zhang  (2004).  which  covers  roughly  the  same  time 
frame.  This  period  ts  long  enough  for  us  to  detect  dec¬ 
adal  changes  in  the  system  by  averaging  over  two  full 
ENSO  cycles  (each  5-6  yr  long).  Choosing  this  period 
to  study  also  allows  us  to  compare  the  model  simulation 
with  high  accuracy  in  situ  ocean  observations,  which  are 
available  for  this  most  recent  decade. 

3.  Results 

a .  STC  pycnoctme  transport  and  its  variability 

The  mean  overturning  circulation  in  the  model's  up¬ 
per  ocean  (Fig.  1)  is  obtained  by  interpolating  meridi- 
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Fra  2.  (a)  Modeled  mean  SSI  averaged  for  all  months  from  1992  to  2003.  (b)  Modeled  SST  difference  between  the  two  periods 
(Period  2  minus  Period  1).  (e)  Observed  mean  SST  averaged  for  alt  months  from  1992  to  2003.  fd)  Observed  SST  difference  between 
the  two  periods  (Period  2  minus  Period  l).  Units  are  °C.  Overplotted  is  the  mean  wind  stress  [in  (a)]  and  the  wind  stress  difference  for 
Period  2  minus  Period  I  [in  (b)]«  Unit  of  wind  stress:  N  m  . 


onal  velocity  tn  HYCOMs  hybrid  vertical  coordinate 
to  fixed  depth  (the  interpolation  interval  is  20  m)  and 
then  integrating  zonal ly  and  vertically  the  interpolated 
meridional  velocity*  The  advantage  of  this  presentation 
is  that  it  gives  a  visual  impression  of  where  the  circula¬ 
tion  takes  place  in  latitude  and  depth:  the  disadvantage 
of  depth  averaging  though  is  that  it  can  create  spurious 
circulation  cells  in  regions  of  sloping  isopycnals 
(Hazeleger  et  al.  2001).  For  example,  the  so-called 
tropical  cells,  which  are  local  circulation  within  5°  of  the 
equator  in  both  hemispheres,  are  a  knowfn  artifact  of 
fixed-depth  averaging.  Ignoring  the  tropical  cells*  circu¬ 
lation  associated  with  subtropical  subduction  is  cen¬ 
tered  between  50-  and  200-m  depths.  This  circulation 
has  intensified  bv  2-3  Sv  (Sv  =  10*  nr  s”1)  in  both  the 
northern  and  southern  branch  from  the  time  between 
July  1992  and  June  1998  (hereinafter  referred  to  as  Pe¬ 
riod  1)  to  between  July  1998  and  June  2003  (hereinafter 
referred  to  as  Period  2)  (Fig.  lc)*  These  two  periods 
w^ere  chosen  as  in  MePhaden  and  Zhang  (2004)  because 
each  spans  one  complete  ENSO  cycle.  In  addition,  the 
observed  Pacific  decadal  oscillation  (PDO)  index 
dropped  significantly  from  Period  I  to  Period  2  (see 
MePhaden  and  Zhang  2004  for  more  details). 


The  simulated  mean  SST  (Fig*  2a)  is  reasonably  close 
to  observations  (Fig.  2c).  In  particular,  the  eastern 
equatorial  cold  tongue  is  captured  in  the  model*  al¬ 
though  its  amplitude  is  somewhat  weaker  than  ob¬ 
served.  Over  most  of  the  Pacific,  the  absolute  values  of 
discrepancies  between  the  simulated  and  observed  SST 
averaged  for  the  full  study  period  are  less  than  I DC  (not 
shown).  The  modeled  SST  difference  between  Periods 
1  and  2  (Fig.  2b)  is  characterized  by  an  antiphase  rela¬ 
tion  between  the  midlatitude  North  Pacific  and  tropical 
Pacific  anomalies.  During  Period  2  when  the  STC  is 
stronger,  positive  SST  anomalies  (relative  to  those  dur¬ 
ing  Period  1 )  developed  in  the  midlatitude  bell  between 
20°  and  40  N,  with  maximum  values  situated  in  the 
middle  of  the  basin:  concomitantly,  negative  SST 
anomalies  exist  in  the  Tropics  and  in  the  near¬ 
shore  waters  off  the  west  coast  of  North  America.  The 
modeled  SST  anomaly  pattern  is  similar  to  its  observa¬ 
tional  counterpart  (Fig.  2d),  Both  of  these  patterns  rep¬ 
resent  the  canonical  pattern  of  SS  I'  anomalies  associ¬ 
ated  with  the  PDO  (Mantua  et  al,  1997:  Lalif  and  Bar¬ 
nett  1994).  although  in  the  model  the  tropical  SST 
anomaly  penetrates  farther  to  the  west  than  in  obser¬ 
vations. 
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Ffg.  3.  Pycnocline  volume  transport  ai  9<;N  and  9*S  as  a  function  of  time.  The  results  at  9°N  are  inverted  for  easy 
comparison,  (a),  (b)  Monthly  data  and  (c)*  (d)  annual  averages  of  the  monthly  means  are  shown.  Black*  red,  and 
green  lines  represent  the  interior  ocean,  western  boundary*  and  total  across  the  basin  transport*  respectively.  The 
interior  ocean  and  western  boundary  subdomain  is  separated  by  145*E  ( IfiO'E)  for  the  transport  at  9°N  (TS). 


The  wind  stress  anomaly  from  Period  1  to  Period  2 
(Fig.  2b)  shows  an  overall  strengthening  of  the  trade 
winds.  Stronger  trade  winds  would  cause  the  tropical 
ther  module  to  slope  down  even  more  to  the  west.  The 
resulting  thermocline  lilt  anomalies  would  be  consistent 
with  the  stronger  equatorward  flow  in  the  pycnocline 
shown  in  Fig,  1c.  Notice  also  that  the  wind  stress 
anomalies  extend  to  the  extratropics  and  may  impact 
subtropical  subduction  and  therefore  tropical  thermo¬ 
cline  ventilation. 

The  volume  transports  in  the  model’s  pycnocline 
across  9°N  and  9°S  as  a  function  of  lime  are  presented 
in  Fig.  3.  They  are  calculated  using  modeled  meridional 
velocity  and  layer  thickness  on  the  original  vertical  co¬ 
ordinates  without  interpolation.  The  latitude  9°N  is 
chosen  because  it  is  a  choke  point  in  the  communica¬ 
tion  between  the  subtropics  and  Tropics  owing  to  the 
existence  of  a  potential  vorticity  (PV)  ridge  at  this  lati¬ 
tude  (Figs.  4a.b);  9°S  is  chosen  for  hemispheric  symme¬ 
try*  Throughout  this  paper*  the  model  pycnocline  is  de¬ 
fined  as  the  layers  between  the  22.5  and  25.3  <xH  surfaces 
(red  lines  in  Fig.  1)  because  ii  captures  the  density 
range  over  which  the  majority  of  equatorward  flow 
takes  place  in  the  model*  This  range  of  density  classes 
significantly  overlaps  with,  but  is  slightly  lower  than. 


that  analyzed  by  McPhaden  and  Zhang  {2002,  2004) 
from  observations*  The  total  transport  across  the  basin 
is  divided  between  the  interior  ocean  (east  of  145°E  at 
9*N  and  east  of  160°E  at  9°S)  and  the  western  boundary 
(longitudes  west  of  145°  and  I  W  E.  respectively).  The 
anticorrelation  between  the  components  (named  as  "in¬ 
terior''  and  "wb"  in  Fig.  3)  suggests  that  transport 
anomalies  through  each  region  tend  to  compensate 
each  other  on  inlerannual  time  scales*  but  there  are 
exceptions*  For  example,  in  the  annual  average  trans¬ 
port  time  series  at  9°S  (Fig.  3d),  the  ant  [correlation 
between  the  western  boundary  and  interior  transport 
deteriorated  in  the  latter  part  of  the  record*  Moreover* 
this  compensation  is  only  partial  such  that  the  trend  in 
the  total  transport  across  the  basin  follows  that  of  the 
interior  ocean  transport.  This  partial  compensation  is 
consistent  with  indirect  estimates  of  western  boundary 
current  (WBC)  transport  (McPhaden  and  Zhang  2004) 
and  other  forced  model  simulations  (Lee  and  Fukumori 
2003;  Hazeleger  et  a!,  2004;  Capolondi  et  aL  2005), 
which  showed  that  the  reduction  in  interior  ocean 
transport  convergence  from  the  mid-1970s  to  earlier 
1990s  was  partially  compensated  by  WBC  transports* 
The  total  pycnocline  transport  convergence  across 
the  basin  averaged  for  Periods  I  and  2  is  1 93)  ±  (Lb  Sv 
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Fr;,  4.  Absolute  value  of  potential  vorticity  (10  lf>  m  1  s  1 )  on  the  24.5-cr„  surface  averaged 
for  (a)  Period  I  and  (b)  Period  2.  (c)  Difference  in  the  mean  depth  of  the  24.5-cr0  surface  for 
Period  2  minus  Period  1. 


and  23.9  ±  0-8  Sv,  respectively.  Standard  errors  through¬ 
out  this  study  are  calculated  as  o/\/7v,  where  a  is  the 
standard  deviation  of  monthly  transport  for  each  period 
and  /V  is  the  number  of  degrees  of  freedom  estimated 
assuming  the  transport  is  an  autoregressive  process 
(Davis  197b).  Therefore*  the  total  pycnocline  transport 
convergence  in  the  model  increased  by  4*9  ±  ID  Sv 
from  Period  1  to  Period  2  (Table  I)*  The  primary  con¬ 
tribution  to  this  increase  comes  from  interior  ocean 


transport  convergence*  which  is  increased  by  4.5  ±  1.6 
Sv  from  3.4  ±  1.2  Sv  in  Period  1  to  7.9  ±  LI  Sv  in 
Period  2  (Table  1).  These  modeled  interior  ocean  mean 
transport  values  are  significantly  lower  than  the  geo- 
strophic  estimates  of  13*4  ±  L6  Sv  in  Period  1  and 
24.1  ±  1.8  Sv  in  Period  2  found  in  McPhaden  and 
Zhang  (2004).  There  are  several  possible  reasons  for 
this  model-observation  discrepancy,  which  will  be  dis¬ 
cussed  later.  However,  most  important  for  our  pur- 
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Table  1 .  Mean  transport  of  various  branches  of  the  STC  and 
the  standard  errors  estimated  from  monthly  data.  Pycnocline 
transports  include  the  transport  convergences  at  9*N  and  9°S  for 
the  interior  ocean  and  total  across-basin  components,  EUC  trans¬ 
port  is  defined  at  140  W.  The  upwelling  is  integrated  over  the 
region  9°N-9°S,  90" W  -180°.  We  estimated  the  standard  errors  of 
the  difference  between  the  means  as  V  e;  +  e%.  where  e,  and  e2  is 
the  standard  error  for  each  period. 


Period  2  - 

Branch 

Period  1 

Period  2 

Period  1 

Total  pycnocline 

19,0  ±  0.6  Sv 

23,9  ±  0.8  Sv 

49  ±  1,0  Sv 

transport 

Interior  pycnocline 

3,4  ±  1,2  Sv 

7,9  ±  LI  Sv 

4,5  £  1.6  Sv 

transport 

EUC 

27.2  ±2.1  Sv 

33.1  ±  2.2  Sv 

5.9  ±  3.0  Sv 

Upwelling 

1 1.7  ±  2.1  Sv 

16.6  i  4.0  Sv 

4.9  ±  4.5  Sv 

poses,  there  is  a  significant  increase  in  interior  ocean 
pycnocline  volume  transport  from  Period  l  to  Period  2, 
in  qualitative  agreement  with  the  observations* 

As  mentioned  earlier,  the  spatial  pattern  of  the  SST 
anomaly  suggests  that,  when  the  STCs  are  stronger 
(weaker),  tropical  SSTs  are  cooler  (warmer)  (Fig.  2b). 
This  relationship  is  further  investigated  in  the  temporal 
domain  (Fig,  5),  The  simultaneous  correlation  coeffi¬ 
cient  between  monthly  interior  ocean  pycnocline  trans¬ 
port  convergence  across  9  N  and  9"$  and  SS  I  anomaly 
averaged  in  the  eastern  tropical  Pacific  (9°N-9°S,  180°- 
90° W)  is  —0.7  (statistically  significant  at  a  95%  confi¬ 
dence  level).  Lagged  cross-correlation  function  (not 
shown)  indicates  that  the  maximum  correlation  occurs 


when  the  former  leads  the  latter  by  2  months.  The  cor¬ 
relation  coefficient  decreases  to  0.4  when  the  total 
transport  across  the  entire  latitude  circle  is  used  in  the 
calculation,  which  is  still  statistically  significant  but  only 
marginally.  This  decrease  in  correlation  is  expected 
given  the  partial  compensation  of  the  western  boundary 
transport  anomalies  and  interior  ocean  transport  varia¬ 
tions. 

The  zonal  structure  of  STC  pycnocline  transport  can 
be  seen  in  ihe  accumulated  integration  of  the  transport 
westward  from  the  eastern  boundary  (Fig.  6).  We  have 
done  this  calculation  for  ihe  9  N  and  9  'S  sections,  as 
well  as  for  the  13DN  and  13°S  sections,  to  investigate  if 
there  are  meridional  variations  in  the  zona)  structure, 
and  also  because  other  models  have  shown  that  trans¬ 
port  variations  at  13°N  and  13  8  are  larger  than  at  9  N 
and  9*S  (Capoiondi  et  ah  2005).  For  the  Northern 
Hemisphere  branch,  most  of  the  equator  ward  flow  oc¬ 
curs  between  140°E  and  181)",  while  for  ihe  Southern 
Hemisphere  branch,  the  equatorward  transport  starts 
farther  east  at  140  W.  These  zonal  structures  are  simi¬ 
lar  qualitatively  for  the  two  latitudes  within  each  hemi¬ 
sphere,  but  the  mean  equatorward  transport  as  well  as 
the  transport  difference  between  Periods  1  and  2  are  in 
general  larger  at  the  higher  latitudes. 

Differences  in  ihe  zonal  structure  between  the 
Northern  and  Southern  Hemispheres  are  caused  by  the 
existence  of  a  potential  vorticity  ridge  in  the  Northern 
Hemisphere  bul  not  in  the  Southern  Hemisphere  (Fig. 
4).  Following  MePhaden  and  Zhang  (2004),  potential 
vorticity  is  defined  as  p“  \fdpidz,  where /is  the  Coriolis 


Fits.  5.  Suite!  line:  Time  series  of  the  interior  ocean  pycnocline  transport  convergence  across 
9£  N  and  ^FS.  Shown  are  monthly  anomalies  relative  to  the  climatological  annual  cycle.  Dashed 
line:  time  series  of  monthly  SST  anomaly  (  C)  averaged  in  the  eastern  tropical  Pacific  from  9  S 
to  9  N  „  I  KM  t  o  90  W ,  Tit  c  SST  time  series  i  s  in  ve  rt  e  d  for  e  asy  c<  j  m  pari  son . 
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. Period  1  — . Period  2  - all 

FiG,  6.  Zonally  accumulated  pyenocline  transport  as  a  function  of  longitude.  The  integration  is  started  from  the 
eastern  boundary  and  accumulates  westward.  The  latitude  at  which  the  computation  is  performed  is  marked  on  the 
plots.  Dashed  lines  represent  averages  for  Period  1,  dotted  lines  for  Period  2.  and  solid  lines  for  all  12  years  of  data. 


parameter,  dpldz  is  the  vertical  density  gradient,  and  pa , 
is  a  constant  reference  density  (10*  kg  m  3).  We  calcu¬ 
late  PV  on  the  24.5-cr»  surface,  which  is  situated  verti¬ 
cally  at  the  center  of  the  model  pyenocline.  The  spatial 
structure  of  the  simulated  PV  field  is  remarkably  close 
to  that  shown  in  McPhaden  and  Zhang  (2004,  their 
Fig.  2),  although  quantitative  differences  still  exist.  The 
spatial  pattern  indicates  that  under  the  influence  of  the 
intertropicai  convergence  zone  (ITCZ),  a  high  PV  ridge 
extends  from  the  Central  American  coast  to  about 
160° W  in  the  northern  subtropics.  Therefore  water  sub¬ 
ducted  in  the  northern  subtropics  has  to  move  west¬ 
ward  parallel  to  PV  contours  (assuming  PV  is  con¬ 
served)  before  turning  southward  in  the  central  basin, 
causing  weak  equatorward  flows  in  the  eastern  basin.  In 
contrast,  water  in  the  Southern  Hemisphere  can  move 
equatorward  even  in  the  eastern  basin  because  of  the 
lack  of  PV  ridge  there. 

The  difference  in  the  mean  depth  of  the  trn  =  24.5 
surface  between  the  two  periods  (Fig.  4c)  shows  that 
the  zonal  slope  of  this  isopycnal  surface  increases  with 
time,  owing  to  shoaling  in  the  eastern  basin  and  deep¬ 
ening  in  the  western  basin  at  9°N,  and  primarily  to 
shoaling  in  the  eastern  basin  at  9°§f  This  increase  in 


isopycnal  slope  is  consistent  with  changes  in  surface 
wind  stress  (Fig.  2b)  and  the  upward  trend  in  the  pyc- 
noeline  transport  (Fig.  3). 

The  separation  of  the  lines  representing  averages  for 
Periods  1  and  2  in  the  interior  basin  and  the  tendency 
for  them  to  merge  in  the  western  boundary  (Fig.  6) 
again  demonstrates  the  partial  compensation  between 
interior  and  western  boundary  transports,  but  the  de¬ 
gree  of  compensation  is  both  temporally  and  spatially 
varying.  For  example,  compensation  is  not  obvious  for 
9°S  (Fig.  6d),  where  the  difference  in  the  mean  trans¬ 
ports  for  the  two  periods  in  the  interior  ocean  remains 
in  the  western  boundary.  Judging  from  the  transport 
time  series  (Fig.  3d),  the  breakdown  of  the  anticorrela¬ 
tion  between  the  interior  and  western  boundary  trans¬ 
port  anomalies  only  occurred  in  the  later  period.  Lee 
and  Fukumori  (2003)  argue  that  the  compensation  is 
caused  by  the  fact  that  the  pyenocline  transport  is  con¬ 
trolled  by  both  the  overturning  circulation  and  the  hori¬ 
zontal  gyre  circulation,  while  Capotondi  el  al.  (2005) 
suggest  that  it  can  be  caused  by  westward  propagation 
of  Rossby  waves  as  the  ocean  dynamically  adjusts  to 
surface  wind  forcing  changes.  Whatever  the  mecha¬ 
nism,  our  results  using  an  eddv-resolving  model  in  gen- 
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pycnocline  _ volume  chcmge  rote 

totol  surface 


Fjg,  7.  Mass  budget  of  ihe  tropical  control  volume  defined  as 
the  region  from  the  surface  to  the  25.3-cre  isopycnal  in  the  vertical 
and  to  9°S  lalitudmally  across  ihe  entire  basin.  Shown  arc 
annually  averaged  data:  ''surface,'"  “pycnocline."  and  “total"  in 
the  legend  refers  to  meridional  transport  convergence  across  the 
basin  through  the  surface  layers,  the  pycnocline  layers,  and  all 
layers,  respectively,  of  the  control  volume. 


eral  corroborate  these  earlier  coarse-resoluiiun  model 
results  about  the  importance  of  western  boundary  cur¬ 
rent  compensation  on  decadal  time  scales. 

h ,  Influence  of  STC  on  the  mass  balance  of  the 
tropical  ocean 

To  examine  the  mass  balance,  we  define  a  control 
volume  as  the  region  from  the  surface  to  the  25.3-<th 
isopycnal  in  the  vertical  and  from  9°S  to  9°N  latitudi¬ 
nally  across  the  entire  basin.  The  total  transport  con¬ 
vergence  across  the  northern  and  southern  boundaries 
of  this  box  is  the  summation  of  the  pycnocline  trans¬ 
port  convergence  and  the  surface  transport  divergence 
(Fig,  7).  As  in  the  above  subsection,  the  pycnocline 
transport  is  calculated  using  model  meridional  velocity 
and  layer  thickness  across  the  basin  between  22.5  and 
25.3  <xB  surfaces  excluding  the  mixed  layer  when  the 
22.5-f ri}  surface  outcrops.  The  surface  transport  diver¬ 
gence  is  calculated  from  the  surface  to  the  22.5-cr^  sur¬ 
face,  or  to  the  base  of  the  mixed  layer  (defined  as  the 
depth  at  which  temperature  is  0.2°C  colder  than  SST) 
when  the  22,5  uH  outcrops.  For  clarity,  we  only  show 
annual  mean  data  in  Fig.  7.  hut  the  conclusion  remains 
the  same  when  monthly  data  is  used.  Both  pycnocline 
convergence  and  surface  divergence  show  a  decadal 
trend  toward  increasing  values  with  time.  These  in¬ 


creases  nearly  compensate  each  other,  causing  the  total 
meridional  iransport  convergence  across  the  lateral 
boundaries  of  the  control  volume  (black  line  in  Fig,  7) 
to  be  quasi  steady  and  much  smaller  in  magnitude  than 
either  component.  As  has  been  observed  in  previous 
modeling  studies  of  El  Nino  and  La  Nina  (e.g..  Springer 
ct  al.  1990),  there  is  also  a  tendency  for  compensation 
between  the  pycnocline  transport  convergence  and  sur¬ 
face  transport  divergence  on  interannual  time  scales, 
especially  between  1997  and  2000, 

Total  meridional  transport  convergence  into  the  con¬ 
trol  volume  is  negative,  and  there  is  a  near-constant 
offset  between  this  convergence  and  the  time  rate  of 
change  of  the  control  volume.  This  offset  implies  a 
steady  upwelling  through  the  base  of  25.3-crw  surface 
between  9°S  and  9°N,  However,  because  this  upwelling 
is  located  at  the  base  of  the  pycnocline*  it  does  not 
represent  a  contribution  of  STCs  to  equatorial  up- 
welling, 

c.  EUC,  equatorial  upwelling,  and  their 

relationships  to  STC  meridional  transport 

In  this  section  we  examine  variations  in  EUC  trans¬ 
port  and  equatorial  upwelling  relative  to  other  branches 
of  the  STCs.  As  a  first  step,  we  validate  the  model  EUC 
against  observations.  The  meridional  and  vertical  struc¬ 
ture  of  the  zonal  flow  in  the  model  (Fig.  8a)  is  very 
similar  to  shipboard  observations  (Johnson  ct  al.  2001. 
their  Fig.  2a).  At  136nW.  the  core  of  EUC  is  located  at 
0°,  120  m  in  the  simulation  (compared  to  0°.  1 10  m  in 
observation),  and  its  maximum  velocity  reaches  90 
cm  s  *,  the  same  as  in  observations.  Results  from  year 
1999  are  shown  here  as  an  example,  but  comparisons 
for  other  years  are  equally  good.  The  simulation  also 
captures  other  observed  current  structures  in  the  tropi¬ 
cal  ocean:  the  northern  branch  of  the  South  Equatorial 
Current  (SEC)  with  a  surface  core  at  2  'N,  the  southern 
branch  of  SEC  with  a  surface  core  at  4°S.  and  the  North 
Equatorial  Countercurrent  (NECC)  with  its  core  cen¬ 
tered  al  50  m  along  7°N.  Although  not  particularly  ob¬ 
vious  in  Fig.  8a,  two  subsurface  countercurrents  also 
exist  in  the  simulation,  at  least  for  some  years.  In  the 
zonal  direction,  the  modeled  EUC  axis  along  the  equa¬ 
tor  tilts  upward  from  west  to  east  (Fig.  8b),  with  the 
maximum  velocities  situated  between  155°  and  120'  W. 
Both  tilt  and  zonal  structure  are  consistent  with  obser¬ 
vations  (Johnson  et  al,  2001,  their  Fig.  3a). 

The  zonal  flow  al  the  equator  has  been  monitored 
continuously  by  ADCP  measurements  from  the  early 
1990s  to  the  present  (Plimpton  el  al.  2004),  providing  an 
Opportunity  for  us  to  compare  these  data  with  the  simu¬ 
lation.  At  14CFW,  vertically  integrated  zonal  flow  in  the 
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Fit.,  8.  Annual  mean  zonal  velocity  field  in  the  model  Results 
are  shown  for  year  1999:  (a)  Meridional  section  at  136*W  from  8°S 
to  ItrN;  (b)  zonal  section  at  the  equator  from  I7(T  to  95°W. 
Contour  interval  is  10  cm  s  1  in  both  plots.  Positive  values  indi¬ 
cate  eastward  flow, 

model  bears  a  striking  resemblance  to  the  A  DCF  mea¬ 
surements  in  both  amplitude  and  temporal  evolution 
(Fig,  9a)*  Comparisons  at  other  longitudes  where  long 
records  of  ADCP  data  are  available  (1  10*W,  170C!W. 
]65°E)  indicate  that  this  good  comparison  is  robust 
across  all  longitudes.  We  then  computed  EUC  trans¬ 
port  in  the  model  as  the  integration  of  zonal  velocity 
with  magnitude  larger  than  10  cm  s  1  from  the  shallow¬ 
est  model  layer  to  400-m  depth  and  from  2DS  to  2°N. 
These  transport  estimates  are  highly  correlated  with 
estimates  of  the  vertically  integrated  zonal  flow  at  the 
equator  (Fig.  9b).  as  in  observations  (c.g,,  Knox  and 
Ha! pern  1982),  Based  on  Figs,  9a  and  9b,  we  conclude 
that  EUC  transports  in  the  model  during  this  period  are 
very  close  to  observations.  The  modeled  EUC  mean 
transport  increased  from  27.2  ±  2,1  Sv  in  Period  1  to 
33.1  ±  2,2  Sv  in  Period  2  {Fig.  9c),  The  derived  differ¬ 
ence  in  the  EUC  transport  (5.9  ±  3.0  Sv)  is  consistent 
with  the  increase  in  the  pycnocline  transport  conver¬ 
gence  across  the  basin  from  Period  I  to  Period  2.  as  well 
as  the  transport  convergence  increase  in  the  interior 
ocean  (see  discussion  in  section  3a  and  Table  1). 

Now  we  compare  STC  pycnocline  transport  conver¬ 
gence  in  the  interior  ocean  with  up  welling  in  the  easl- 


50 lid:  ADCP.  dashed:  HYCOU 


20,  60,  TOO.  )  40.  180  220. 


Fig.  9.  (a)  Zonal  flow  at  thu  equator  integrated  from  the  shal¬ 
lowest  data  point  to  400-m  depth  (units:  nr  s  solid  line  from 
ADCP  data  and  dashed  from  the  HYCQM  simulation.  The  shal¬ 
lowest  (deepest)  ADCP  data  at  this  longitude  are  at  KJ-m  (415-m) 
depth,  (b)  Scatterplot  of  modeled  zonal  flow  ai  the  equator  vs  the 
modeled  EUC  transport  (Sv).  Zonal  flow  integrated  both  verti¬ 
cally  and  meridionally  (from  2  S  to  2"N).  The  vertical  integration 
in  both  cases  is  from  the  shallowest  model  layer  to  400-m  depth, 
(cl  EUC  transport  time  series  from  the  model:  Horizontal  lines 
denote  the  mean  transport  for  Period  1  and  Period  2,  respectively. 
Ail  data  shown  are  at  140  \V\ 


ern  tropical  Pacific.  Using  the  total  across-basin  trans¬ 
port  did  not  change  the  basic  conclusion  although  it 
lowered  the  correlation  coefficient  (see  details  below). 
Here  the  upwelling  transport  is  calculated  as  the  spa¬ 
tially  integrated  vertical  velocity  at  the  base  of  the 
mixed  layer.  The  spatial  domain  is  chosen  to  be  9f,N- 
9°S,  180  -90  "W  [the  same  index  region  for  SST  vari¬ 
ability  in  McPhaden  and  Zhang  (2004)  and  Fig.  4 
herein J,  The  resulting  upwelling  transport  is  highly  cor¬ 
related  with  the  pycnocline  transport  convergence  with 
no  obvious  lag  (Fig.  10).  The  simultaneous  correlation 
coefficient  between  the  monthly  time  series  (Fig.  10a)  is 
0.76  (statistically  significant  at  a  95%  confidence  level, 
lowered  to  (164  when  cross- basin  pycnocline  transport 
convergence  was  used).  The  averaged  upwelling  trans¬ 
ports  for  Periods  1  and  2  are  11.7  ±  2.1  Sv  and  16.6  ±  4.0 
Sv,  respectively.  These  values  correspond  well  writh  the 
1993-99  mean  cross-isothermal  upwelling  transport  near 
the  base  of  the  surface  mixed  layer  of  15  Sv  computed 
from  observed  winds  and  upper-ocean  thermal  field  be¬ 
tween  5qN  and  5qS,  95°  and  155nW  (Meinen  et  al,  2001 ), 
The  interrelationship  between  pycnocline  transport 
in  the  interior  ocean,  EUC,  and  eastern  basin  upwelling 
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Fig.  UK  (a)  Relationship  between  the  interior-ocean  pycnodine 
transport  convergence  across  9° N  and  9°S  (solid  line,  using  the  left 
coordinate:  units:  Sv)  and  the  upwelling  at  the  base  of  mixed  layer 
spatially  integrated  over  9DN-9 'S.  1 80  -90  W  (dashed  line,  using 
the  right  coordinate;  units:  Sv).  (b)  As  in  (a),  except  now  the 
monthly  data  were  filtered  by  a  13-point  running  mean  filter. 

is  summarized  in  Fig,  !  I.  These  time  series  are  highly 
cross  correlated  at  zero  lag.  The  correlation  coefficient 
between  the  filtered  pycnodine  transport  convergence 
and  EUC  transport  (upwelfing)  is  0,7  (0.6).  which  is 
statistically  significant  with  a  95%  confidence  level 

The  absolute  mean  transport  amplitudes  of  these 
variables,  in  addition  to  those  of  the  total  transport 
convergence  across  the  basin  for  Periods  1  and  2*  and 
Lhe  differences  between  the  two  periods,  are  listed  in 
Table  1.  Bearing  in  mind  lhe  uncertainties  represented 
by  the  standard  errors  in  each  transport  estimate,  these 
values  indicate  that  the  increase  in  upwelling  between 
the  two  periods  (4.9  ±  4.5  Sv)  is  dynamically  consistent 
with  the  increases  in  both  equator  ward  pycnodine 
transport  and  EUC  transport;  furthermore,  the  increase 
in  the  across- basin  pycnodine  transport  convergence  is 
dominated  by  that  in  the  interior  ocean. 

d.  Impacts  on  meridional  heat  transport 

To  understand  the  impact  of  STC  variations  on  the 
upper-ocean  heat  balance,  we  calculated  meridional 
heat  transport  associated  with  STC  as 


H(y,t)=  ^ 

2=pyc  J 


<  VhT)\IyM dx 


-  TEk(y 


.o2  f 

Z  [>VC  J 


(  vh  Av.y.z.!  dx.  (1) 


Fick  1L  Comparison  between  interior-ocean  pycnodine  trans¬ 
port  convergence  (black).  EUC  transport  (solid  line),  and  the 
eastern  basin  upwelling  magnitude  (United  line).  For  each  vari¬ 
able,  a  monthly  climatology  based  on  all  12  years  of  data  for  each 
variable  was  removed:  the  resulting  monthly  anomaly  time  series 
was  then  subjected  to  a  13-poinl  running  mean  filter.  Last*  lhe 
filtered  data  were  normalized  by  their  respective  standard  devia¬ 
tions,  The  corresponding  standard  deviation  for  the  pycnodine 
transport  convergence,  EUC  transport,  and  upwelling  magnitude 
is  3,8,  6.4,  and  4,1  Sv.  respectively. 


where  E^=pyc/  dx  denotes  zonal  integration  across  the 
basin  and  then  vertical  summation  over  the  range  of  the 
pycnodine,  and  w  h.  and  T  are  layer  meridional  veloc¬ 
ity,  thickness,  and  temperature,  respectively.  Here  // 
can  be  scaled  to  PW  by  multiplying  with  pCp  x  It) 
where  p  =  1026  kg  m  1  and  Cp  =  399(1  J  K  1  kg 

Equation  (I)  represents  an  idealized  model  of  heat 
transport  bv  the  STCs  in  which  we  assume  that  the 
upper  branch  of  the  STC  in  the  Ekman  layer  carries  the 
same  amount  of  transport  as  in  the  pycnodine  layer 
with  a  zonal  mean  temperature  of  TEk(  v*  /)♦  In  practice* 
rEk(y,  t)  is  calculated  as  the  volume-weighted  zonal- 
mean  mixed  layer  temperature.  The  advantage  of  this 
formula  is  that  it  ensures  a  mass  balance  so  that  we  can 
define  a  “heal  transport"  associated  with  the  STCs  as  a 
function  of  time.  Errors  introduced  by  using  this  for¬ 
mula  will  be  discussed  later. 

Thus  defined,  the  long-term  mean  heat  transport  di¬ 
vergence  of  //  (computed  by  subtracting  heat  transport 
at  9CS  from  that  at  9°N)  is  0.67  ±  0,04  PW,  in  good 
agreement  writh  the  heat  transport  divergence  associ¬ 
ated  with  the  overturning  circulation  using  another  ide¬ 
alized  model  and  National  Centers  for  Environmental 
Prediction  (NCEP)  surface  wind  data  (Hazeleger  el  al. 
2004), 

Because  vh  and  T can  be  written  as  vft  ~  (vh)  +  (vh)1 
and  T  =  T  +  T\  where  (  )  represents  long-term  lime 
averages  and  a  prime  represents  deviations  from  the 
long-term  means,  the  total  heat  transport  //  can  be  de¬ 
composed  as 
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F*o.  12.  Relationship  between  pycnodine  volume  transport 
convergence  across  9'  N  and  9°S  (solid  lines,  units:  SvT  using  the 
right  >  axis)  and  heat  transport  divergence  across  the  same  lati¬ 
tudes  (dashed  lines  lines,  units:  PW.  using  the  left  y  axis,  negative 
means  convergence).  Heal  transport  divergence  is  decomposed 
into  contributions  from  anomalous  circulation  and  thermal 
anomalies  (named  as  HI.  H2,  and  H3)  according  to  Eqs.  (2)-(4) 
in  the  text.  The  correlation  coefficient  between  the  two  time  series 
in  each  panel  is  marked. 

H  -  HJvh*  T.  7\Tk)  +  H}(vh\  T,  7^) 

+  H2(vtK  T\  rBk)  +  H3(vh\  T\  T'Bk)x  (2) 

where  //,„  represents  the  time-invariant  part  of  //.  //, 
represents  the  effect  of  circulation  anomalies  adverting 
the  mean  thermal  field,  //:  represents  the  effect  of 
mean  circulation  advecting  thermal  anomalies,  and  H3 
represents  the  effect  of  circulation  anomalies  advecting 
thermal  anomalies. 

The  divergence  of  Ht,  H2 ,  and  //3  across  9°N  and  TS 
and  its  relationship  to  pycnodine  transport  conver¬ 
gence  across  the  same  latitudes  is  shown  in  Fig.  12.  The 
scale  of  the  left-hand  y  axis  in  Fig.  12a  is  approximately 
one  order  of  magnitude  larger  than  that  in  Fig.  12b. 
Hence,  variability  of  the  total  heat  transport  divergence 
is  predominantly  controlled  by  /Y,.  which,  not  surpris¬ 
ingly,  is  highly  correlated  with  pycnodine  transport 
convergence  anomalies  at  zero  lag  (r  =  0.77,  which  is 
statistically  significant  at  95%  confidence  level).  In 
other  words,  volume  transport  variability  explains 
about  60%  (r  —  0,59)  of  the  variance  in  heat  transport. 
Interestingly.  ff2  is  anticorrelated  with  pycnodine 


!0G*E  HCPE  180°  140°W  10OQW 


Fig.  13.  Net  surface  hem  flux  (W  m  ’).  positive  means  into  the 
ocean,  (a)  averaged  between  1992  and  2003.  and  (b)  the  difference 
Period  I  minus  Period  2.  Contour  interval  is  III  (5)  for  (a)  [(b)]. 


transport  convergence  anomalies  (r  -  —0.42,  which  is 
marginally  significant  at  95%  confidence  level).  The 
negative  correlation  can  be  explained  by  the  interaction 
of  ocean  circulation  and  property  fields:  as  the  circula¬ 
tion  intensifies  from  Period  I  to  Period  2,  tropical  up¬ 
per-ocean  temperature  decreases,  causing  the  heat 
transport  divergence  associated  with  the  thermal 
anomalies  to  decrease:  hence,  a  downward  trend  of  H2 
is  seen  in  Fig.  12b,  The  magnitude  of  the  eddy  compo¬ 
nent  H3  is  even  smaller  than  H2 
The  notion  that  the  upper-ocean  thermal  balance  is 
perturbed  by  lateral  heat  transport  convergence/ 
divergence  is  supported  by  concomitant  changes  in  the 
net  surface  heat  flux  (Fig.  13).  From  Period  I  to  Period 
2,  net  surface  heat  flux  into  the  ocean  between  9°S  and 
9°N  has  increased  by  0,08  ±  0.05  PW  (equivalent  to 
2.3  ±  1.4  W  m  UC  '),  while  the  Ht  heat  transport 
divergence  across  9°S  and  9°N  has  increased  by  0.13  ± 
0,04  PW.  The  cooling  of  the  upper  ocean  from  Period  1 
to  Period  2  is  therefore  caused  by  the  dominance  of 
anomalous  lateral  heat  transport  divergence  (a  cooling 
effect)  over  net  surface  heat  flux  anomalies  (a  heating 
effect).  Using  observations  and  different  OGCM  simu¬ 
lations,  Liu  and  Huang  (2000)  attributed  the  1945-93 
SST  warming  trend  in  the  tropical  Pacific  to  ocean 
advection.  In  their  case,  tropical  SST  increases  while 
net  surface  heal  flux  has  a  cooling  trend  of  roughly 
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10  W  m  :  (50  vr)  Both  the  negative  feedback  be¬ 
tween  SST  and  net  surface  heal  flux  and  the  trend  of 
net  surface  heat  flux  [2,3  ±  1.4  W  m  ?  ( 12  yr)  1  in  our 
study  versus  10  W  m  ?  (50  yr)  1  in  Liu  and  Huang 
(2000)]  are  consistent  beiween  the  two  studies. 

The  0.13  PW  increase  in  the  heat  transport  diver¬ 
gence  of  Hx  is  20%  of  the  long-term  mean  of  the  total 
heat  transport  divergence  associated  with  STC,  which  is 
0.67  PW.  It  is  also  comparable  to  the  anomalous  heat 
transport  divergence  caused  by  extra  tropica  I  wind 
anomalies  in  an  OGCM  coupled  to  a  mixed-layer  at¬ 
mospheric  model  and  to  El  Nino-induced  anomalous 
overturning  heat  transport  divergence  in  that  same 
model  (Haze lege r  el  al.  2004).  In  other  words,  this  mag¬ 
nitude  of  change  in  the  tropical  ocean  heal  transport 
divergence  seems  robust  across  different  models. 

The  main  sources  of  error  in  the  idealized  formula  of 
Eq.  (1 )  include  the  use  of  zonal  mean  temperature  for 
calculating  temperature  transport  in  the  Ekman  layer. 
Also,  in  our  model  the  contribution  of  thermal  anoma¬ 
lies  to  the  heat  transport  variations  is  likely  under¬ 
estimated  because  spiciness  anomalies  associated  with 
surface  freshwater  forcing  are  suppressed.  Because  of 
these  errors,  the  exact  values  of  //  and  its  components 
should  be  viewed  with  caution,  but  the  trends  are  ro¬ 
bust  and  can  be  explained  by  the  underlying  physics.  As 
the  STC  increases,  the  associated  poleward  heal  trans¬ 
port  increases,  causing  the  heal  transport  out  of  the 
Tropics  to  increase;  hence  the  cooling  of  the  tropical 
upper  ocean  occurred  despite  increasing  net  surface 
heat  flux  into  the  ocean  at  the  same  time.  As  the  tropi¬ 
cal  ocean  cools,  the  poleward  heat  transport  associated 
with  temperature  anomalies  decreases,  but  the  overall 
cooling  is  determined  by  the  intensifying  circulation. 
The  result  that  circulation  anomalies  dominate  over 
temperature  anomalies  in  bringing  about  poleward  heal 
transport  variations  is  also  found  in  a  100()-yr  integra¬ 
tion  of  a  global  coupled  model  (Merry field  and  Boer 
2005). 

4.  Discussion 

a.  Pycnodine  transport  and  ext  rat  topical  processes 

Since  pycnodine  transport  in  our  calculations  is  de¬ 
fined  at  relatively  tow  latitudes  (9°N  and  9nS),  it  is  pos¬ 
sible  that  this  transport  and  its  variations  are  only 
driven  by  tropical  winds.  However,  we  cannot  rule  out 
the  involvement  of  higher  latitudes  based  on  a  two- 
dimensional  Eulerian  description  of  the  circulation.  To 
examine  how  the  simulated  pycnodine  transport  in  the 
Tropics  is  related  to  extralropical  subduclion.  we  cal¬ 
culated  the  subduclion  rate  as  the  horizontally  inte¬ 
grated  vertical  velocity  at  the  base  of  the  seasonal  ther- 


mocline  defined  by  the  maximum  winter  mixed  layer 
depth.  We  neglected  contributions  from  lateral  induc¬ 
tion,  which  were  found  to  be  small  for  this  region 
(Huang  and  Qiu  1994,  1998).  The  winter  (January- 
March  for  the  North  Pacific  and  June-August  for  the 
South  Pacific)  subduclion  rale  averaged  over  the  1 2 -yr 
studying  period  is  approximately  8.5  Sv  (nr  37  m  yr  '  ) 
in  the  South  Pacific  between  10°  and  20  S.  160  and 
100  W  and  1 1.8  Sv  (or  45  m  yr  1 )  in  the  North  Pacific 
between  10  and  35  N.  170  and  140  W.  The  difference 
of  the  subduclion  rate  between  Periods  !  and  2  (Period 
2  minus  Period  I )  ts  3.0  Sv  and  0.9  Sv  for  the  North  and 
South  Pacific,  respectively.  These  numbers  are  consis¬ 
tent  with  the  modeled  pycnodine  transport  variations 
shown  in  Table  1.  With  the  approximation  used  in  our 
calculation,  they  are  also  broadly  consistent  with  obser¬ 
vations. 

To  ultimately  demonstrate  the  influence  of  extralro¬ 
pical  processes  on  the  tropical  ocean  circulation,  one 
needs  to  use  Lagrangian  methods.  Using  a  passive 
tracer  simulation.  Fukumori  el  al.  (2004)  showed  that 
roughly  70%  of  the  water  mass  of  the  Niho-3  region  can 
be  traced  back  to  eastern  subtropical  ihermodine  wa¬ 
ters  on  decadal  time  scales:  furthermore,  eddy  transport 
due  to  intraannual  variability  significantly  enhances  the 
transport  magnitude  of  the  interior  ocean  pathways. 
Passive  tracer  simulations  are  beyond  the  scope  of  this 
study,  but  they  would  be  crucial  for  understanding  the 
three-dimensional  pathways  connecting  the  extra  tropi¬ 
cal  and  tropical  ocean  circulations, 

b ,  Relationship  between  pycnodine  transport 
and  SST 

A  study  by  Capotondi  el  al.  (2005)  using  the  National 
Center  for  Atmospheric  Research  (NCAR)  OGCM 
suggests  that  SST  anomalies  in  the  equatorial  up  welling 
region  lead  the  STC  transport  convergence  anomaly 
across  9°N  and  9  S  by  2  months  at  maximum  correla¬ 
tion  (their  Fig.  14).  whereas  in  this  model  SST  anoma¬ 
lies  lag  the  STC  transport  convergence  anomalies  by 
2  months  (Fig.  14).  In  both  models  the  correlation 
remains  statistically  significant  all  the  way  across  the 
basin  interior,  indicating  that  pycnodine  transport 
changes  arc  relevant  for  SS T  changes  in  those  regions. 
The  difference  in  sign  of  the  corresponding  phase  rela¬ 
tionship  between  the  two  models  needs  to  be  studied 
further  because  the  role  of  STC  transport  in  regulating 
tropical  SST  depends  on  the  relative  timing  of  each 
process.  Capotondi  et  al  attributed  the  time  lag  of 
pycnodine  transport  relative  to  SST  anomalies  in  the 
NCAR  model  to  zonal  averaging  of  transport  values. 
Because  the  STC  transport  adjusts  from  the  east  to  the 
west  via  oceanic  Rosshy  waves,  zonal  averaging  across 


May  2007 


CHENG  ET  At. 


1353 


1 7Q®£  1 70aW  1 50°W  MG 1  IOBW 


Flo.  14.  Cross-correlation  between  zonal  I  y  accumulated  (starl¬ 
ing  at  KH)  W  and  ending  at  160  E)  interior-ocean  pycnocluie 
transport  convergence  anomalies  across  9  S  and  9°N  and  SST 
anomalies  averaged  for  9  N  9°S,  1 80  -90°  W.  Positive  (negative) 
lags  mean  SST  anomaly  lags  (leads)  the  transport  convergence 
anomaly.  Shaded  areas  indicate  that  the  correlations  are  statisti¬ 
cally  significant  at  a  95%  confidence  level  with  the  number  or 
degrees  of  freedom  estimated  assuming  autoregression  process 
(Davis  1976).  Dashed  (solid)  lines  indicate  negative  (positive)  cor¬ 
relation.  Contour  interval  is  (1.1. 

the  basin  can  make  the  STC  transport  appear  to  lag  (be 
SST  anomaly,  which  is  averaged  over  the  central  and 
eastern  Pacific, 

We  repeated  the  lagged  correlation  between  the  zon¬ 
al  ly  accumulated  pycnocline  transport  ami  eastern 
equatorial  SST  anomaly  as  in  Capolondi  el  al.  (2005) 
using  our  model  data.  The  tendency  of  the  maximum 
correlation  to  shift  in  the  direction  of  greater  SST  lead 
relative  to  the  pycnocline  transport  variation  as  one 
moves  westward  is  also  present  in  our  model  (Fig.  14) 
between  170  and  130'  W,  but  in  the  western  basin  (west 
of  170  W)  this  effect  is  not  as  obvious  as  that  in  the 
NCAR  GGCM.  The  overall  smaller  zonal  slope  of  the 
maximum  correlation  axis  in  this  model  compared  to 
that  in  the  NCAR  OGCM  reflects  the  model  depen¬ 
dence  of  this  result.  One  possible  reason  for  such  model 
sensitivity  is  that  the  baroclinic  Rossby  wave  speed  is 
different  between  the  two  models.  Although  a  detailed 
comparison  of  baroclinic  Rossby  wave  speed  between 
the  models  is  beyond  the  scope  of  this  study,  we  specu¬ 
late  that  the  HYCOM  baroclinic  wave  speed  is  gener¬ 
ally  larger  (because  its  across- basin  adjustment  time  is 
quicker  than  the  NCAR  model  according  to  Fig.  14) 
owing  to  its  smaller  vertical  diffusivity  in  the  pycnocline 


and  stronger  stratification.  Such  model  sensitivities 
need  to  be  kept  in  mind  when  interpreting  results  of 
model  simulations. 

Nevertheless,  both  model  simulations  show  that, 
when  only  transport  in  the  eastern  part  of  the  basin  is 
considered,  the  pynocline  transport  convergence  varia¬ 
tion  leads  the  SST  changes  by  a  few  months.  There  are 
two  possible  explanations  for  this:  1)  eastern  tropical 
SST  is  in  part  responding  to  perturbations  associated 
with  pycnocline  transport  anomalies;  2)  both  eastern 
tropical  SST  and  pycnocline  transport  are  responding 
to  the  same  forcing,  such  as  surface  wind  anomalies,  but 
the  pycnocline  transport  responds  quicker  than  SST, 
The  exact  reason (s)  for  this  lead-lag  relationship  are 
still  not  precisely  understood,  hut  our  results  indicate  a 
close  relationship  between  SSI1  and  transport  variations 
on  interannual  to  decadal  lime  scales. 

c.  Model  limitations 

The  model  outputs  analyzed  in  this  study  are 
monthly  averages,  so  the  effects  of  submonthly  time 
scale  variability  Is  not  accounted  for.  High  frequency 
variability  may  induce  additional  mass  transport  varia¬ 
tions  due  to  correlations  between  velocity  and  isopyc- 
nal  layer  thickness  on  time  scales  shorter  than  1  month. 
However,  for  the  domain  studied  here,  the  biggest  eddy 
effects  are  restricted  to  the  equatorial  band  within  5  of 
the  equator  (Hazeleger  et  al.  2001).  Therefore,  we  do 
not  expect  our  results  to  be  significantly  different  when 
submonthly  scale  eddy  effects  are  ineluded. 

We  also  have  to  keep  in  mind  the  limited  model  ge¬ 
ometry,  Giese  el  al.  (2002)  found  a  strong  influence  of 
South  Hemisphere  processes  on  the  tropical  and  North 
Pacific  thermodine  temperature  anomalies.  Our  model 
uses  a  buffer  zone  al  the  southern  boundary  to  simulate 
the  South  Pacific  inflow  to  the  domain.  Although  it  has 
been  suggested  that  the  overturning  circulation  is  rela¬ 
tively  insensitive  to  the  buffer  layer  setup  al  the  lateral 
boundaries  in  an  isopycnal  model  (e,g„  Smith  et  al. 
2000),  the  variability  of  pycnocline  transport  across  the 
southern  boundary  in  our  model  is  likely  underesti¬ 
mated,  as  reflected  tn  the  smaller  than  observed  geo- 
strophic  transport  differences  between  Periods  1  and  2 
(section  3a),  Also,  the  Indonesian  Through  flow  (1TF) 
passage  is  closed  in  the  model.  However,  correspon¬ 
dence  of  our  results  with  those  of  Lee  and  Fukumori's 
(2003)  suggests  that  our  basic  conclusions  may  he  in¬ 
sensitive  to  whether  the  1TF  is  present  or  not.  Never¬ 
theless.  Hazeleger  et  al.  (2004)  found  that  the  1TF  is 
important  for  the  heat  transport  compensation  between 
the  overturning  circulation  and  gyre  circulation  in  their 
model:  more  specifically,  when  the  ITF  is  reduced, 
poleward  heal  transport  in  the  South  Pacific  weakens. 
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The  precise  role  of  the  ITF  in  balancing  the  mean  mass 
and  heat  transport  between  different  basins  on  decadal 
time  scales  needs  to  be  further  clarified. 

In  addition  to  model  geometry,  potential  biases  in 
model  physics  could  affect  the  comparison  with  obser¬ 
vations  as  well.  For  example,  McCreary  and  Lu  ( 1994) 
and  Lu  et  aL  (1998)  found  that  the  PV  ridge  in  a  layer 
model  tends  to  be  too  strong  and  causes  an  interior 
ocean  equatorward  flow  that  is  too  weak.  In  addition, 
differences  can  arise  because  of  the  sensitivity  of  mod¬ 
els  to  the  specification  of  mixing  para  me  terizat  ions. 
Uncertainties  in  ocean  properties  specified  at  lateral 
open  boundaries  and  in  surface  forcing  variables  and 
computational  errors  all  may  have  contributed  to  dis¬ 
crepancies  between  transports  in  the  model  and  obser¬ 
vations.  For  instance,  surface  wind  fields  differ  greatly 
among  different  products  in  the  tropical  Pacific  (e.g., 
Hackerl  et  aL  2001)  and  these  differences  can  lead  to 
large  differences  in  patterns  and  strength  of  oceanic 
circulation  (McPhaden  and  Zhang  21K12,  2004).  Uncer¬ 
tainties  in  surface  buoyancy  fluxes  are  also  very  large  in 
the  tropical  Pacific  and  the  specification  of  these  fluxes 
can  affect  model  results  (Wang  and  McPhaden  2001). 
In  view  of  these  limitations,  the  model  has  performed 
reasonably  well  in  capturing  the  correct  sense  of  the 
increased  intensity  of  the  circulation  in  the  tropical  Pa¬ 
cific  from  Period  I  to  Period  2. 

5.  Conclusions 

The  main  conclusions  from  this  study  are  summa¬ 
rized  below: 

*  Consistent  with  observations,  the  HYCOM  simula¬ 
tion  shows  a  trend  toward  increasing  pycnoclinc  vol¬ 
ume  transport  across  the  basin  from  1992  to  2003. 
This  increase  is  attributed  primarily  to  changes  in  the 
interior  ocean  pycnodine  transports  while  variability 
at  the  western  boundary  partially  compensates  that  in 
the  interior.  However,  this  partial  compensation  is 
not  operative  at  all  times  and  at  all  latitudes. 

*  Meridional  volume  transport  convergence  anomalies 
in  the  model  subtropical  pycnodine  are  in  quantita¬ 
tive  agreement  with  anomalous  volume  transports  in 
both  the  observed  and  modeled  Equatorial  Undercur¬ 
rent,  as  well  as  modeled  anomalous  upwelling  into  the 
mixed  layer  in  the  eastern  tropical  Pacific  (Table  1). 

*  Consistent  with  the  increasing  circulation  intensity, 
lateral  heat  transport  out  of  the  tropical  control  vol¬ 
ume  (defined  as  the  region  between  9°N  and  9US,  and 
from  the  surface  to  it,,  =  25.3  isopycnal)  increases 
from  Period  i  to  Period  2.  leading  to  a  cooling  of  the 
tropical  upper  ocean  despite  that  net  surface  heat 


flux  into  the  control  volume  has  increased  over  the 

same  lime. 

The  overall  purpose  of  this  study  is  to  shed  light  on 
ocean  physics  that  control  tropical  Pacific  responses  to 
imposed  surface  forcing  variability  using  an  eddy- 
resolving  general  circulation  model.  We  found  that  the 
decadal  changes  in  pycnodine  transport  convergence, 
EUC.  equatorial  upwelling,  and  SST  in  the  eastern 
tropical  Pacific  vary  in  concert  with  one  another.  As 
such,  changes  in  STCs  and  the  associated  lateral  heat 
transport  divergence  anomalies  play  a  central  role  in 
regulating  the  tropical  Pacific  climate  on  decadal  time 
scales. 

The  ocean-only  model  used  in  l his  study  does  not 
allow  it  to  be  used  to  understand  ocean-atmosphere 
interactions,  and  the  eddy-resolving  resolution  prohib¬ 
its  sensitivity  experiments,  in  particular,  we  have  not 
addressed  how  equatorial  SST  anomalies  may  influence 
the  atmosphere  circulation  and.  furthermore,  how  the 
atmospheric  feedbacks  may  affect  STC  transports. 
Coarser- resolution  ocean  model  experiments  suggest 
that  both  tropical  and  extratropical  wind  anomalies  can 
have  an  impact  on  the  STCs  and  tropical  SST  variability 
on  interannual  and  decadal  time  scales  (Nonaka  el  aL 
2002).  These  forcing  mechanisms  are  likely  to  play  a 
role  even  in  eddy-resolving  models. 

The  model  and  the  observations  are  in  qualitative 
agreement  on  several  key  aspects  of  oceanic  variability. 
But.  there  are  still  significant  quantitative  differences 
between  the  modeled  interior  pycnodine  transport  val¬ 
ues  and  observations,  as  discussed  earlier.  Apparently, 
our  ability  to  better  understand  and  more  accurately 
simulate  the  earth's  climate  and  its  variations  relies  on 
further  improvement  in  the  design  and  implementation 
of  global  ocean  and  atmosphere  observation  systems, 
belter  estimates  of  air-sea  fluxes,  and  advances  in 
model  physics  formulations. 
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